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ABSTRACT 


This  report  summarizes  the  work  done  by  the  SDL  during  the 
period  July  through  September  1968,  and  is  primarily  concer  led 
with  seismic  research  activities  related  to  the  detection  and 
identification  of  nuclear  explosions  and  earthquake  phenomena. 
Also  discussed  are  the  support  tasks  and  data  services  per¬ 
formed  for  other  participants  in  the  VELA-Uniform  projects. 
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I.  INTRODUCTION 

This  quarterly  report  summarizes  the  technical  work,  sup¬ 
port  effort,  and  service  tasks  completed  during  the  period 
July  through  September  1968.  Current  or  past  work  is  mentioned 
only  if  it  relates  to  the  present  discussions. 

Reviews  of  technical  reports  completed  during  the  report¬ 
ing  period  are  contained  in  Section  II  under  descriptive 
headings.  Section  III  is  a  summary  of  the  support  and  service 
tasks  performed  for  in-house  projects  and  for  other  VELA- 
TTniform  participants.  This  report  concludes  with  Appendix  A 
in  which  are  listed  the  organizations  who  received  SDL  data 
services  during  the  period. 

H .  WORK  COMPLETED 

A .  Preliminary  Evaluation  of  NORSAR 

The  purpose  of  this  study  was  to  present  a  preliminary 
evaluation  of  tne  Norwegian  Seismic  Array  (NORSAR) .  We  were 
interested,  specifically,  in  determining:  1)  the  minimum  inter- 
sensor  spacing  required  to  produce  optimum  rms  noise  reduction 
by  summing  the  array,  and  2)  the  amount  of  signal  loss,  rms  noise 
reduction,  and  signal-to-noise  gain  produced  by  beamforming  the 
array.  The  basic  procedures  included  reformating  the  original 
data,  prefiltering,  computing  power  spectra,  time-shifting,  and 
summing. 


The  results  of  the  analysis  of  one  noise  sample  are 
included  here  and  are  summarized  in  Figure  1.  The  rms  noise 
reduction  produced  by  summing  the  various  combinations  of  sub- 
arrays  is  shown  as  a  function  of  the  corresponding  average  inter- 
sensor  spacing  and  compared  with  the  corresponding  N^  (shown  by 
the  dashed  lines) .  These  results  show  that  to  obtain  optimum 
rms  noise  reduction,  relative  to  N^,  by  prefiltering  (0.7  -  5.0 
cps)  and  summing  the  array,  the  minimum  inter-sensor  spacing 
must  be  approximately  three  kilometers. 

Figures  in  the  original  report  (No.  221)  show  the 
results  of  the  analysis  of  the  noise  power  spectra  for  the 
various  combinations  of  subarrays  used  to  form  the  unphased 
sums.  In  these  figures,  the  power  spectra  of  the  unphased  sum 
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RMS  NOISE  REDUCTION  (db) 
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AVERAGE  INTER-SENSOR  SPACING  (km) 

Figure  1.  RMS  noise  reduction  by  summing  various  combinations 
of  NORSAR  elements 
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of  the  noise  sample  was  compared  with  the  average  of  the  spectra 
from  all  of  the  subarrays  used  to  form  that  sum.  These  figures 
show  that  the  spectral  shape  is  essentially  the  same  for  both 
the  unphased  sum  and  the  average  of  the  individual  spectra.  The 
figures  also  show  that,  with  the  exception  of  the  example  where 
the  average  inter-sensor  spacing  is  only  2.5  km,  the  unphased 
sum  attenuates  the  noise  at  1  cps,  for  example,  by  a  factor 
comparable  to  N*5.  The  spectral  analyses  were  performed  in  all 
cases  on  data  prefiltered  to  the  band  0.7  -  5.0  cps. 

Table  1  summarizes  the  results  of  the  amplitude  analysis 
of  the  two  teleseismic  events  under  three  main  headings:  Signal, 
Noise,  and  S/N.  The  amount  of  signal  loss,  rms  noise  reduction, 
and  signal-to-noise  gain  produced  by  beamforming  the  array  using 
all  12  elements  which  were  operational  at  the  time  the  events 
were  recorded  is  shown  under  the  corresponding  heading.  As 
shown  in  Table  1,  the  signal  loss  produced  by  beamforming  the 
array  is  less  than  one  db  for  both  events.  The  rms  noise  re¬ 
duction  achieved  by  forming  the  beams  is  equal  to  the  N*5 
criterion  for  both  events.  The  signal-to-noise  gain  from  the 
beamsteered  array  is  attenuated  one  db  relative  to  for  the 
two  events. 

The  preliminary  evaluation  of  NORSAR  short-period  data 
included  an  analysis  of  one  signal-free  noise  sample  and  the 
records  from  two  teleseismic  events.  The  following  conclusions 
were  based  on  the  results  of  that  study: 

1.  The  minimum  inter-sensor  spacing  required  to  produce 
an  optimum  rms  noise  reduction,  relative  to  N^,  by  prefiltering 
and  summing  the  array  outputs  is  approximately  three  kilometers. 

2.  A  signal  loss  of  less  than  one  db  is  produced  in  beam- 
forming  the  array  by  time-shifting  the  subarray  data  using  only 
computed  travel-time  differences.  Part  of  this  signal  loss 
results  from  using  the  Fc  filter.  The  remainder  of  the  signal 
loss  is  attributed  to  small  misalignment  of  the  P  wave  and  to 
differences  in  wave  form  across  the  array. 

3.  In  all  cases  where  the  average  inter-sensor  spacing  is 
greater  than  three  kilometers,  the  rms  noise  reduction  is  equiv¬ 
alent  to  the  criterion.  This  same  noise  reduction  factor  is 
shown  by  the  power  spectra  at  1  cps. 
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4.  A  3-to-l  ratio  of  Range-to-RMS  is  obtained  in  all  cases 
for  a  100-second  noise  sample. 

5.  The  signal- to -noise  ratio  produced  by  the  beamformed 
array  is  attenuated  one  db  relative  to  N^5. 

Reference 


Clark,  D.M. ,  1968,  "Preliminary  Beamforming  Study  of  the  TFO-37 
Array",  Report  No.  216,  Seismic  Data  Laboratory,  Teledyne,  Inc. 
Alexandria,  Virginia. 

B*  A  Preliminary  Evaluation  of  the  Matched  Filter  Technique 
iQ _ the  Detection  of  Long— Period  Body  Wave  Radiation. 

SDL  Report  No.  222  concerns  the  detection  of  long- 
period  body  wave  radiation  for  both  synthetic  and  real  test  cases 
using  the  matched  filter  approach. 

To  determine  the  effectiveness  of  the  matched  filter 
technique,  we  first  set  up  synthetic  test  cases  by  burying  a 
known  earthquake  signal  in  long-period  noise  at  various  S/N 
ratios.  We  then  used  the  known  earthquake  signal  as  a  matched 
filter.  An  optimum  filter  length  was  determined  from  the  results 
of  this  experiment,  along  with  an  estimate  of  the  minimum  S/N 
ratio  at  which  the  technique  breaks  down.  The  results  from  the 
synthetic  cases  were  sufficiently  encouraging  to  warrant  proceeding 
the  search  for  long-period  body  waves  from  actual  events.  For 
actual  events  the  results  are  less  encouraginy.  We  attribute  this 
m  part  to  the  lack  of  sufficient  available  data  from  LASA.  Only 
LASA  data  was  considered  because  it  is  already  in  digital  format 
and  is  best  suited  for  application  of  array  summing  techniques. 

The  first  sections  of  the  report  deal  with  the  methods 
and  results  for  the  synthetic  test  cases,  and  the  latter  sections 
eal  with  the  application  of  the  technique  to  observed  data. 

The  results  obtained  in  this  study  are  mixed.  The 
synthetic  test  cases  show  that  when  the  signal  used  as  the  fil¬ 
ter  is  identical  with  the  unknown  signal, the  matched  filter 
technique  alone  is  successful  in  detecting  the  unknown  signal 
at  a  minimum  S/N  ratio  of  4.  By  beamforming  the  matched  filter 
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outputs  from  the  nine  long-period  sensors  in  the  A,  E,  and  F  rings 
at  LASA,  we  were  able  to  detect  the  presence  of  the  unknown  signal 
at  S/N  ratios  as  low  as  1.  The  much  poorer  performance  of  the 
matched  filter  on  observed  data  may  be  attributed  to  differences 
within  the  wavetrain  of  the  reference  event  with  respect  to  the 
unknown  event,  and/or  the  relative  absence  of  long-period  body 
wave  radiation  from  the  smaller  events. 

The  matched  filter  technique  worked  best  on  the  seismo¬ 
grams  from  a  Kurile  Island  double  event.  Using  the  whole  suite 
of  arrivals  P  through  S,  we  were  able  to  detect  both  events,  the 
second  with  an  apparent  polarity  reversal.  However,  the  S/N 
ratio  of  the  matched  filter  sum  trace  for  the  first  event  was 
considerably  less  than  that  obtained  by  straightforward  beam¬ 
forming  the  band  passed  seismograms.  This  suggests  that  there  is 
considerable  mismatch  in  waveform  between  the  reference  signal 
and  the  unknown  signals,  even  though  the  events  occurred  close 
together  within  the  same  source  region.  Also,  the  marked  improve¬ 
ment  in  S/N  ratio  using  the  shorter  500-second  filter  suggests 
that  the  mismatch  becomes  more  important  when  the  later  phases, 
particularly  S  arrivals,  are  included  in  the  filter.  The  second 
event  is  not  easily  detected  from  the  individual  traces,  but 
after  beamforming  the  matched  filter  outputs, it  is. 

Three  events  recorded  on  April  1,  1967,  had  epicenters 
within  55  km  of  one  another.  Furthermore,  the  epicenters  reported 
by  the  USC&GS  are  all  at  the  same  depth,  the  P  wave  delays  across 
LASA  being  almost  identical.  Therefore,  the  observed  signal  mis¬ 
match  is  probably  due  to  differences  in  relative  excitation  of  the 
various  arrivals  between  the  events.  Theoretically,  we  would 
expect  this  difference  to  be  greatest  for  the  S  wave  arrivals. 
However,  we  were  not  able  to  eliminate  the  effect  of  the  differ¬ 
ences  using  the  level  or  rectified  filters,  which  suggest  that 
the  differences  may  be  more  complex  than  simple  variations  in 
amplitude  or  polarity  reversals  of  the  arrivals. 

The  matched  filter  results  from  the  December  22,  1966 
event  were  interesting.  The  focal  depth  for  this  event  was  77  km. 
It  is  located  some  350  km  from  the  reference  event,  although  the 
P  wave  delays  across  LASA  for  the  two  events  are  equal  to  within 
0.5  seconds.  Using  the  730-second  filter,  which  includes  the  P 
thr  agh  SKS  arrivals,  we  were  able  to  detect  lcng— period  body  wave 
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radiation  at  the  appropriate  times  for  both  P  and  PcP.  This  is 
probably  due  to  the  fact  that  the  S-P  interval  is  of  very  nearly 
the  same  duration  as  the  SKS-PcP  interval  at  this  particular 
distance  and  focal  depth.  The  S/N  ratio  of  the  matched  filter 
outputs  is  low,  and  without  prior  knowledge  of  the  expected 
arrival  time  it  is  doubtful  whether  the  onset  would  have  been 
recognized,  particularly  using  the  500-second  filter. 

In  the  case  of  three  events  with  magnitudes  less  than 
5,  we  were  unable  to  detect  any  long-period  body  wave  radiation, 
either  by  matched  filtering  or  beamforming.  This  in  part  may  be 
due  to  differences  in  epicenter  location,  particularly  as  they 
a^e  at  different  depths  from  the  reference  events.  However,  the 
partia1  success  of  the  matched  filter  on  the  77  km  deep  event 
suggests  that  more  probably  these  events  radiated  only  small 
amounts  of  energy  in  the  pass  band  of  10-50  seconds. 

The  failure  of  the  matched  filter  below  magnitude  5 
leads  to  the  question  of  threshold.  Obviously  the  main  problems 
involved  are  radiation  pattern  differences  and  differences  in 
polarity  from  one  event  to  another.  Relative  time  shifts  of 
arrivals  between  the  reference  event  and  the  unknown  event  will 
also  tend  to  degrade  the  matched  filter  performance,  since  the 
wavetrains  no  longer  match  exactly.  Differences  in  radiation 
pattern  result  not  only  in  a  mismatch  of  the  wavetrains,  but 
also  in  differences  in  reported  body  wave  magnitudes.  Through¬ 
out  the  report  we  gave  magnitudes  based  on  the  USC&GS  Bulletin, 
the  North  American  stations  reported  in  this  Bulletin,  and  the 
magnitudes  reported  from  LASA.  The  variation  indicates  the 
difficulty  in  assigning  a  threshold.  Our  results  (Figure  2) 
indicate  that  the  threshold  for  detection  is  somewhere  around 
mb  =  5  for  the  Kurile  Islands  in  terms  of  USC&GS  magnitudes. 

It  is  impossible  to  express  this  result  in  terms  of  LASA  magni¬ 
tudes  since  the  LASA  magnitude  (m£>)  of  .ne  reference  event  as 
reported  is  less  than  5. 

Finally,  we  used  the  portions  of  the  filter  between 
the  principle  arrivals  as  our  filter  in  several  instances. 

These  intervals  are  the  ones  in  which  PL  modes  are  most  likely 
to  occur.  We  found  no  evidence  that  these  sections  alone  im¬ 
proved  our  ability  to  detect  long-period  energy.  Nor  did  simply 
shortening  this  interval  to  bring  the  principle  phases  in  the 
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Log  S/N  Ratio  (J 9  Channels) 


USC&GS  Short  Period  Body  Wave  Magnitude 


Figure  2.  Graph  showing  S/N  ratios  obtained  by  matched 

filtering  as  a  function  of  body  wave  magnitude 

(mb) 


filter  into  alignment  with  the  predicted  arrival  times  of  the 
unknown  event  produce  any  significant  improvement  in  S/N  ratio 
for  the  smaller  (m^  <  5)  events. 

The  long-period  body  wave  matched  filter  technique 
although  it  is  not  as  effective  as  the  surface  wave  technique' 
m  terms  of  detecting  whether  or  not  an  event  occurred  in  a 
particular  region,  may  constitute  an  important  tool  for  the 
classification  of  events.  If  any  long-period  body  waves  are 
found,  their  presence  would  indicate  that  the  event  is  an 
earthquake,  since  long-period  teleseismic  P  waves  are  seldom 
observed  for  explosions.  Moreover,  the  method  is  cheap  in 
terms  of  computer  time  and  could  be  incorporated  in  an  on-line 
detection  system.  The  ultimate  goal  should  be  a  system  in 
which,  once  an  event  has  been  detected  and  located  by  other 
means,  a  long-period  body  wave  matched  filter  could  be  quickly 
and  inexpensively  applied  to  determine  the  presence  or  absence 
of  long-period  body  wave  radiation. 

From  the  results  of  this  study  we  conclude: 

1.  For  cases  where  the  unknown  signal  and  the  reference 
signal  are  very  similar  in  waveform,  the  matched  filter  works 


2-  For  input  S/N  ratios  less  than  4,  array  summing  of  the 
individual  match  filter  outputs  is  required  to  observe  lonq- 
penod  body  wave  radiation. 


c3,;  I°ealli'.  beamforming  the  matched  filter  outputs  irnprov. 
the  S/N  ratio  by  a  factor  of  2  (6  db)  over  beamforming  the  band- 
pass  filtered  seismograms. 

4.  Using  observed  data,  the  threshold  for  detection  of 
long  period  body  waves  from  events  in  the  Kurile  Islands  at 
depths  .ess  than  80  km  is  probably  around  magnitude  5  (idk) 
using  9  elements  of  LASA.  The  best  results  are  likely  to  be 

obtained  when  the  depth  of  the  reference  event  and  the  unknown 
event  are  the  same. 


5.  The  matched  filter 
wave  arrivals  are  used. 


technique  works  best  when  only  the 
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6.  It  would  not  appear  legitimate  to  include  portions  of 
the  seismogram  after  the  S  arrival  time  and  attribute  the  output 
to  the  presence  of  long-period  body  waves,  since  we  know  that 
higher  mode  surface  waves  have  S  wave  velocities  as  an  upper 
bound  and  may  be  present  in  the  record  anywhere  after  the  S  arrival 

7.  Despite  several  attempts  to  use  only  "interphase"  (PL) 
portions  as  a  reference  filter  we  found  no  evidence  that  using 
these  PL  portions  of  the  record  alone  improved  our  ability  to 
discern  the  presence  of  long-period  energy. 

8.  Amplitude  equalization  is  not  an  effective  means  of  im¬ 
proving  the  performance  of  the  matched  filter,  probably  because 
weaker  portions  of  the  reference  signal  contain  noise  which  we 
build  up  by  this  procedure. 

9.  Many  more  events  must  be  analyzed  before  a  reliable 
threshold  for  detection  of  long-period  body  waves  at  LASA  can  be 
established. 
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Principles  of  Wiener  Auto-Adaptive  Filter. 


SDL  Report  N.? .  224  discusses  various  techniques  of 
designing  auto-adaptiv*  (self-changing)  filters,  and  the  appli¬ 
cation  of  one  particular  technique  to  two  teleseisms.  Also 
described  are  the  calculation  of  running  correlation  functions, 

Js mg  a  moving  time  window,  and  iterative  solutions  to  the  multi¬ 
channel  normal  equations.  Moreover  a  listing  is  included  in  the 
original  report  of  a  program  which  will  calculate  an  auto-adaptive 

mu  t2  ~channel  filter  by  the  use  of  running  correlation  functions 
and  the  method  of  steepest  descent. 


Section  III  of  the  report  discusses  the  theory  of  the 
mean  square  error  criterion  and  the  application  of  the  principle 
of  orthoganality  to  multichannel  filtering  as  well  as  various 

adaptive  filter  solutions  which  depend  upcn  running  correlation 
functions. 


A  f il ter  whose  coefficients  depend  upon  the  time- 
varying  statistics  of  the  noise  is  termed  an  adaptive  filter. 

A  filter  which,  as  the  noise  statistics  change,  automatically 
updates  itself  without  the  need  of  human  intervention  or  super¬ 
vision  is  termed  an  auto-adaptive  (self-changing)  filter.  The 
auto  adaptive  filters  used  to  produce  the  following  results 
were  designed  to  be  linear  but  time-varying  multichannel  filters 


Figure  3  is  the  result  of  applying  the  auto-adaptive 
filter  to  the  LASA  input  data  recorded  in  subarray  D1  The 
first  trace  is  a  timing  trace  with  1-second  pips,  every  fifth 
one  of  which  is  accentuated.  Traces  2  through  11  are  the  ob- 
served  data  trace  11  being  the  one  to  be  predicted;  this  trace 
is  shifted  1  second  to  the  left.  Trace  12  is  the  output  of  the 
auto-adaptive  filter.  Trace  13  (the  error  trace)  is  the  dif¬ 
ference  between  trace  11  and  trace  12.  Trace  14  is  the  sum  of 

th!  n^ne  lnPut  traces  used  to  predict  trace  11.  The  output  of 
the  adaptive  filter  is  seen  to  be  identical  to  the  sum  at  the 

,  ^ .  °f  record  (as  lt  should  be,  since  the  initial  guess 

for  the  filters  was  1/N,  0,0,0...). 
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Figure  4  is  a  continutation  in  time  of  the  same  data. 
We  can  see  that  the  filter  is  indeed  learning  from  the  input 
data,  since  the  noise  power  is  getting  progressively  less  in 
trace  12  as  compared  to  trace  14.  In  general  the  same  peaks 
seem  to  be  coming  through,  but  the  power  is  quite  a  bit  less  in 
the  output  of  the  adaptive  filter  (trace  12)  as  compared  to  the 
beamed  sum  trace  (trace  14)  as  time  goes  on.  Figure  5  is  a 
continuation  of  the  same  data  which  includes  the  Aleutian  event. 
The  first  motion  of  the  signal  seems  to  be  enhanced  as  does  a 
secondary  phase,  probably  pP.  The  signal— to-noise  ratio  improve¬ 
ment  is  2  db  over  the  bandpass  phased  sum  due  to  some  signal  re¬ 
jection  by  the  adaptive  filter.  The  second  phase  occurs  roughly 
nine  seconds  after  the  primary  phase,  the  time  delay  being  about 
right  for  pP. 

It  should  be  noted  that  the  progressive  decrease  in 
the  noise  power  with  time  is  due  to  two  factors:  one  is  the 
convergence  of  the  adaptive  filter  from  the  initial  guess  toward 
the  optimum  time-stationary  filter;  the  other  factor  is  the  im¬ 
provement  achieved  by  reducing  time-varying  coherent  noise.  The 
small  improvement  (2  db)  above  the  band  limited  phased  sum  seems 
to  indicate  that  the  amount  of  time-varying  coherent  noise  power 
is  low.  Of  course,  another  factor  to  consider  is  a  poor  array 
response  which  could  cause  coherent  noise  to  alias  back  into  the 
main  lobe. 


From  this  report  we  concluded: 

1.  The  method  of  updating  the  correlation  functions  used  in 
this  report  is  a  valid  method  for  calculating  relatively  stable 
and  slowly  time  varying  correlation  functions. 

2.  The  c'’'' to— adaptive  filters  designed  from  the  time— varying 
noise  correlation  functions  do  converge,  and  they  yield  better 
results  (2  db  better  on  the  LASA  event)  than  simple  bandpass 
filtering  and  beamforming.  The  concept  of  auto-adaptive  filter¬ 
ing  could  be  applied  to  single  channel  (prediction,  prediction- 

error,  and  deghost)  filtering  as  well  as  to  multichannel  filter¬ 
ing. 
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D.  Body  Wave  Magnitude  and  Source  Mechanism 

In  this  study  a  new  method  is  proposed  for  improv¬ 
ing  the  body-wave  magnitude  determination  by  using  the  observed 
values  of  the  body-wave  magnitude  (mb)  together  with  the  first 
motion  directions,  to  obtain  by  least  squares  analysis  the  best 
double  couple  source  parameters;  the  resulting  radiation  pattern 
is  then  integrated  spatially  to  provide  a  corrected  estimate  of 
the  magnitude. 

We  suggest  that  a  better  measure  of  the  "true"  body- 
wave  magnitude  may  be  obtained  by  taking  the  amplitude  of  a 
purely  compressional  source  whose  radiation  pattern  has  the 
same  area  as  the  observed  radiation  pattern. 

Naturally,  this  means  that  we  have  to  know  the  source 
parameters,  i.e.,  the  dip  direction,  dip  angle  and  slip  angle  of 
the  fault  plane  and  the  auxiliary. 

To  find  the  source  parameters  we  propose  to  use  data 
which  are  readily  available,  i.e,  reported  first  motion  directions 
together  with  observed  body-wave  magnitudes  as  given  by  the  C&GS 
Earthquake  Data  Reports  and  the  ISC  Bulletin.  Of  course,  the 
method  is  bound  to  be  inaccurate  if  there  ate  large  gaps  in  the 
azimuthal  coverage,  but  in  this  case  the  other  methods  also  fail. 
It  is  hoped  that  for  earthquakes  of  magnitude  5  and  larger,  an 
approximate  solution  of  the  source  parameters  may  be  obtained  and 
hence  a  better  value  of  the  mean  body— wave  magnitude  may  be  es¬ 
tablished.  For  a  detailed  mathematical  description  of  the  method, 
the  reader  is  referred  to  SDL  Report  No.  225. 

Results  for  a  number  of  events  previously  studied  by 
other  investigators  are  presented  in  the  original  report.  For 
our  purpose  here,  however,  the  analysis  of  one  event,  the  Banda 
Sea  Earthquake  (21  March  1964),  is  described. 

For  this  earthquake  37  stations  were  used,  12  of  them 
with  magnitude  and  25  with  only  first  motions.  Figure  6  shows 
the  observed  radiation  pattern  and  Table  2  shows  the  input  data. 
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Figure  6.  Observed  radiation  pattern  for  the  Banda  Sea 
earthquake  of  21  March  1964 
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Table  2.  Example  of  input  data  for  the  Banda  Sea  earthquake 
of  21  March  1964. 


Figure  7  shows  the  calculated  radiation  pattern  and  Table  3 
gives  the  results  in  tabular  form.  In  our  final  solution, 

36  stations  had  the  same  calculated  and  observed  first  motion 
direction  and  the  one  disagreement  was  a  near  nodal  value. 

The  new  mean  value  of  mb  came  out  to  be  5.85  as  against  the 
C&GS  value  of  5.8.  The  stars  in  Figure  7  are  calculated  am¬ 
plitudes  for  those  stations  which  had  observed  values  and  the 
dots  are  calculated  amplitudes  for  those  stations  which  only 
reported  first  motion  directions. 

In  conclusion,  the  comparison  of  the  fault  plane 
solutions  obtained  by  the  new  method  and  those  obtained  using 
first  motions  only  shows  that  it  is  possible  to  find  reasonable 
agreement.  The  solution  of  the  Hindu-Kush  earthquake  parameters 
show  that  although  better  agreement  with  the  first  motion  di¬ 
rections  was  obtained  by  Hedayati  and  Hirasawa  (1966) ,  the  am¬ 
plitude  radiation  pattern  does  not  agree  as  well  with  the 
observed  pattern  as  that  obtained  by  using  a  combination  of 
first  motion  directions  together  with  observed  body-wave  mag¬ 
nitudes. 


The  results  of  the  Rat  Island  earthquake  of  the  5th 
of  February  1965  show  that  the  proposed  method  compensates 
correctly  where  too  many  of  the  observed  magnitudes  were  near 
to  nodes.  This  shows  that  the  proposed  definition  of  magnitude 
is  superior  to  taking  the  arithmetic  mean  of  the  observations  of 
mb- 


The  fault  plane  solution  found  by  this  method  is 
probably  not  as  accurate  as  that  obtained  by  S-wave  data. 
However,  the  new  method  gives  two  checks  on  the  accuracy 
obtained,  the  first  by  simply  the  number  of  differences  in 
sign  between  calculated  and  observed,  i.e.,  the  "score".  The 
second  is  by  examining  the  confidence  limits  of  the  magnitude. 

The  results  also  indicate  that  in  order  to  take  the 
radiation  pattern  into  account  when  finding  the  "true"  value 
of  mb,  it  is  sufficient  to  use  an  approximate  solution  to  the 
source  parameters.  This  is  possible  by  using  readily  available 
magnitude  and  first  motion  data  without  laboriously  having  to 
re-examine  records. 
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Table  3.  Results  for  the  Banda  Sea  earthquake  of 
21  March  1964. 
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E.  Earthquake  and  Explosion  Analysis. 

Projects  completed  include  photo  albums  for  four 
events  and  travel-time  computations  for  NTS  and  other  sites.  In 
addition  complexities  were  computed  for  a  number  of  special  events. 

III.  SUPPORT  AND  SERVICE  TASKS 


A.  VELA-Uniform  Data  Services 


As  part  of  the  contract  work  statement,  the  SDL  pro¬ 
vided  one  or  more  of  the  following  support  and  service  functions 
for  VSC  and  other  VELA  participants: 

-  copies  of  16  and  35  mm  film 

-  playouts  of  earthquakes  and  special  events 

-  copies  of  composite  analog  tapes 

-  composite  analog  tapes  of  special  events 

-  use  of  1604  computer  for  checking  out  new  programs  or 
running  production  programs 

-  copies  of  digital  programs 

-  digitized  data  in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604. 

-  running  SDL  production  programs,  such  as  power  spectral 
density  and  array  processing  on  specified  data. 

-  digital  x-y  plots  of  power  spectra  or  digitized  data 

-  signal  reproduction  booklets 

-  space  for  visiting  scientists  utilizing  SDL  facilities  to 
study  data  and  exchange  information  with  SDL  personnel. 
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During  this  report  period,  51  such  projects  were 
completed  and  the  10  organizations  receiving  these  services 
are  listed  in  Appendix  A. 

B.  Data  Library 

The  data  library  contains  approximately  8345  digitized 
seismograms,  246  digital  computer  programs  and  306  composite 

analog  magnetic  tapes,  all  available  for  use  by  the  VELA-Uniform 
program. 


The  following  additions  were  made  during  this  period: 
!•  Digital  Seismograms  -  215  including 

-  data  from  three  explosions 

-  eight  earthquakes  recorded  at  various  stations 

2.  LASA  Data  -  30  tapes 

-  there  are  a  total  of  1438  digital  tapes  in  the  library 
including  1034  field  tapes.  There  is  also  a  master 
calibration  tape  which  contains  the  magnification  (digital 
counts  per  millimicron)  of  each  sensor  for  every  subarray. 
These  magnifications  have  been  computed  for  all  calibra¬ 
tion  tapes  currently  in  house.  As  each  new  calibration 

is  received,  it  is  routinely  run  through  the  new  program 
CALIBR  and  added  to  the  master  tape. 

3.  Digital  Programs  -  9  including 

MAPR I NT  -  This  program  works  with  output  produced  by  the  NETWORK 3 
program  on  a  BCD  magnetic  t=»pe.  The  NETWORK3  produces  up  to  9 
different  probability  or  magnitude  values  for  each  event  of  the 
given  set  of  events.  The  purpose  of  MAPRINT  is  to  print  these 
values  so  that  a  specific  map  may  be  overlayed  and  the  values 
contoured.  Up  to  9  different  maps  may  be  printed  corresponding 
to  the  nine  different  values  per  event.  Each  map  contains  all 
events  and  they  have  the  same  overlay. 
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ARRAYTAB  -  Generates  tape  reel  114  which  contains  the  coordinates 
of  arrav  sensors  relative  to  the  center  of  the  array. 

SUB2LIB  -  Uses  all  channels  on  a  subset  tape  to  generate  an  un¬ 
packed  library  format  tape  and  a  listing  cf  the  label  information. 

DIFFER  -  Despikes  seismograms  written  in  the  SDL  SUBSET  format. 

WHITE  -  Finds  a  filter  from  a  selected  time  interval  of  the  first 
trace  of  a  seismogram  and  uses  it  to  filter  all  traces  of  that 
seismogram  over  another  selected  time  interval. 

HRVKSPEC  -  Computes  the  high-resolution  wavenumber  spectrum  of 
seismic  linear-array  time  series  data.  The  high-resolution 
wavenumber  spectrum  may  be  calculated  with  respect  to  a  particular 
sensor  in  the  array,  or  it  may  be  averaged  over  all  of  the  sensors 
in  the  array. 

GEODIST  -  Given  a  reference  point,  GEODIST  computes  the  distance 
(in  both  km  and  degs) ,  azimuth,  and  back  azimuth  to  each  point 
in  an  array  of  secondary  points  using  Rudoe ' s  method  supplied 
by  Mr.  Paul  Thomas  of  NRL. 

NETMEV( modified)  -  This  is  a  modification  of  NETW0RK2  which  is 
designed  to  produce  an  output  tape  containing  event  locations 
and  selected  stations  to  be  processed  by  the  program  LOCMEV. 

LSTSQELP  -  Finds  the  least  squares  fit  to  a  set  of  input  points 
and  computes  a  confidence  ellipse  on  the  slope  and  intercept  of 
the  solution. 

4.  Analog  Composite  Tapes  -  7  including 

a.  Made  by  SDL 

SCROLL 
FAULTLESS 
KNICKERBOCKER 
PIN  STRIPE 
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b.  Made  by  Geo tech 

SHUFFLE 

BOXCAR 

CHATEAIGAY 

C.  Data  Compression 

This  is  a  continuing  routine  operation,  and  production 
is  maintained  at  the  level  needed  to  meet  the  requirements  of 
the  field  operation  (LRSM  and  U.S.  Observatories)  and  the  Seismic 
Data  Laboratory.  For  this  period  1728  tapes  were  compressed. 

D.  Automated  Bulletin  Process 

The  July  1968  Observatory  Bulletin  was  processed 
during  this  report  period  and  forwards’  to  Geotech,  a  Teledyne 
Company,  for  checking  and  publication. 
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APPENDIX  A 


July  -  September  1968 


Organizations  Receiving  SDL  Data  Services 


Massachusetts  Institute  of  Technology 
University  of  Washington 
Texas  Instruments,  Inc. 

St.  Louis  University 

Missile  Electronic  Warfare  Technical  Area,  New  Mexico 
UCLA 

Pennsylvania  State  University 

Air  Force  Office  of  Scientific  Research 

Vitro 

General  Atronics  Corporation 
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